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METHOD AND SYSTEM FOR FILTERING A SIGNAL AND 
PROVIDING ECHO CANCELLATION USING AN 
ADAPTIVE LENGTH FILTER 

5 FIELD OF THE INVENTION 

The present invention relates generally to an adaptive filter with improved operational 
characteristics, and, more specifically, to methods and systems of achieving echo cancellation 
employing the improved adaptive filter. 

10 BACKGROUND OF THE INVENTION 

Adaptive filters are used in numerous applications to remove undesired frequencies from 

O a signal. In an exemplary application, adaptive filters are used in telecommunication systems, 

O 

bj more specifically in echo cancellation systems, to remove from a signal echoes that may arise as 
a result of the reflection and/or retransmission of modified input signals back to the originator of 

fU 15 the input signals. Commonly, echoes occur when signals that were emitted from a loudspeaker 

O 

B are then received and retransmitted through a microphone (acoustic echo) or when reflections of 

jjjj a far end signal are generated in the course of transmission along hybrids wires (line echo). 

Although undesirable, echo is tolerable in a telephone system, provided that the time 
n delay in the echo path is relatively short. But, longer echo delays can be distracting or confusing 

^20 to a far end speaker. Understandably, the telecommunications industry has devoted substantial 
resources to developing systems that minimize echo, without adversely affecting the ability of 
two speakers to communicate with one another. 

Acoustic echo often occurs in speakerphones that employ one or more microphones 
together with one or more speakers to enable "hands-free" telephone communication. Line echo 
25 originates because telephone facilities usually comprise two-wire circuits within each area 

connecting individual subscribers with the switching office and four-wire transmission circuits 
between switching offices in different local exchange areas. A call between subscribers in 
different exchange areas is carried over a two-wire circuit in each of the areas and a four-wire 
circuit between the areas, with conversion of speech energy between the two and four-wire 
30 circuits being conducted by hybrid circuits. If the hybrid circuit input ports had perfectly 
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matched impedances of the two and four- wire circuits, the signals transmitted from one exchange 
area to the other would not be reflected or returned to the first area as echo. Unfortunately, due 
to impedance differences that inherently exist between different two and four-wire circuits, and 
because impedances must be matched at each frequency in the voice band, it is difficult for a 
5 given hybrid circuit to perfectly match the impedances of any particular two and four- wire 
transmission circuit. 

To substantially remove echoes from a communication system, echo cancellation systems 
and methods deploy adaptive filters to generate an estimate (echo estimate) of the echo signal 
that is then removed from the signal being transmitted back to the originator of the echo- 
10 generating signal (far-end source). More specifically, a far-end source transmits a signal (far-end 
signal) that passes through a connection medium and into an input terminal of a communication 
unit. The far-end signal received at the input terminal is cross-coupled via a cross coupling path 
(either acoustically or in line) and creates a cross coupling echo component. That echo 
component combined with a new signal from the near-end (near-end signal) is transmitted back 



m 15 to the far-end source as a composite output signal. In this system, an echo canceller is 



E 



conventionally deployed to monitor the far-end signal and generate an estimate of the actual 
M echoes expected to return in the form of a composite signal with the near-end signal. The echo 

estimates are then applied to a subtracter circuit in the transmit channel to remove or at least 
reduce the actual echo. 

20 To create an accurate estimate of the actual echoes, various types of adaptation methods 

are known in the prior art and can be deployed in the echo canceller in the form of an adaptive 
filter. Conventionally, a finite impulse response (FIR) filter is used that carries out the 
convolution between the far-end signal and the estimated impulse response, N samples in length, 
of the echo paths. In the most basic model, an adaptive filter (filter vector) operates on a far-end 
25 signal vector to produce an estimate of the echo, which is subtracted from the combined near-end 
and echo signal. The overall output of the adaptive echo canceller is then used to control 
adjustments made to tap values of the filter vector. 

In the aforementioned application, and other applications requiring the use of adaptive 
filters, a critical design requirement is the ability for the filter system to achieve convergence in a 
30 rapid, stable manner and, in the process, use a minimal amount of computational resources. A 
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trade-off traditionally exists between stability, accuracy, and speed of convergence of an 
adaptive filter. As such, different types of adaptive filter systems have been deployed in echo 
cancellation systems that attempt to improve the rate of convergence and stability, while still 
minimizing the computational resources required. United States Patent No. 5,995,620 discloses 
5 a method of canceling an echo that, according to the inventors, has an improved convergence 
time with low complexity. The echo cancellation method includes the step of canceling the echo 
in a far-end signal with a Kalman filter having a time varying Kalman gain vector K(t) 
proportional to the vector: 

(p 1 (t)x(t-l)...p n (t)x(t-n)) r 
10 where pi(t) are the diagonal elements of a diagonal matrix P(t) satisfying a Riccati equation, i 
denotes the ith diagonal elements of P(t), t denotes discrete time, n denotes the number of filter 
O taps, and T denotes transpose. 

n 

y Despite the aforementioned prior art, an adaptive filter is still needed that achieves stable 

£8 convergence using fewer computational resources. Additionally, a method and system of echo 

m 

m 15 cancellation displaying stable convergence, while minimizing the computational resources 



required, is also needed. 



H* SUMMARY OF THE INVENTION 

LjjjS 

In 



The present invention is directed toward a novel adaptive filter and novel methods and 
systems for conducting echo cancellation in telecommunication systems. In one embodiment of 
H= 20 the adaptive filter of the present invention, a filter is chosen, preferably an adaptive finite 

impulse response (FIR) filter of an appropriate length N. Once the filter is chosen, convergence 
is achieved using a convergence process. With convergence complete, the total energy of the 
filter coefficients is computed and two energy thresholds, a first energy threshold and a second 
energy threshold, are set relative to the computed total energy. Using the set energy thresholds, 
25 new filter coefficients are determined by relating the starting filter coefficient with the first 
energy threshold and the end filter coefficient with the second energy threshold. A first filter 
delay, in accordance with the determined start coefficient, is incorporated into the filter 
operation. 

In a second embodiment of the present invention, a filter is chosen, preferably an 
30 adaptive finite impulse response (FIR) filter of an appropriate length N. Once the filter is 
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chosen, convergence is achieved using a convergence process, preferably a least means square 
(LMS) approach. With convergence complete, a set of coefficients (Selected Coefficients), 
which represents a sub-set of the entire set of converged coefficients and is preferably pre- 
designated but, optionally, may be dynamically chosen, is used to calculate the average energy of 
5 the portion of the signal corresponding to the Selected Coefficients. The average energy is then 
scaled using a scaling factor to produce a threshold value that is used to select the coefficient 
subset used to filter the signal. 

The present invention is also directed toward novel echo cancellation methods and 
systems that use the novel adaptive filters. Adaptive filters are used to generate an echo signal 
1 0 component used to cancel the echo generated by the engagement of a far-end signal with a cross- 
coupling path. The novel adaptive filter method and system presented herein can be used to 
p improve the calculation of the echo impulse response by, among other things, reducing the 

computational complexity and memory requirements of the coefficient calculation conducted 
S3 within the adaptive filter. In one embodiment, a filter delay is incorporated to insure that only 

1 5 coefficients approximate to or greater than the start coefficient (SC) are employed. In a preferred 
embodiment, the delay value is a function of the SC minus 1 . One of ordinary skill in the art 
would appreciate that different delay values, as related to SC could be used and depends upon on 
how much adherence one of ordinary skill in the art would want the designed system to have 
^ relative to the calculated SC. The modeled response is produced as a function of the transfer 

M= 20 function denoted by z . The values of D2 extend from the point after the initial delay up to a 
point proximate to the end coefficient (EC). In a case where the delay is equal to the start 
coefficient minus 1 and the termination point is chosen to equal the EC, keeping in mind the 
termination point could be slightly greater than or less than EC, depending upon design 
considerations, the transfer function z" can be modeled as: 

25 Z -SC Z -(SC+1) z -(SC + 2) z -(SC + 3)_ z -(EC-3) z <EC-2) Z (EC-I) Z ~EC 

Accordingly, the present invention provides for adaptive filters that have improved 
convergence, computational, and memory bandwidth proprieties. When applied to 
telecommunication applications, the present invention additionally provides for improved 
methods and systems of canceling echoes. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

These and other features and advantages of the present invention will be appreciated as 
they become better understood by reference to the following Detailed Description when 
considered in connection with the accompanying drawings, wherein: 
5 Fig. 1 is a flowchart describing the operation of one embodiment of a novel adaptive 

filter method; 

Fig. 2 is a block diagram of one embodiment of the novel adaptive filter system ; 
Fig. 3 is a flowchart describing the operation of a second embodiment of a novel adaptive 
filter method; 

10 Fig. 4 is a block diagram of a telecommunication system having a voice over packet 

gateway; 

Fig. 5 is a block diagram of an echo cancellation system; and 

Fig. 6 is a block diagram of a novel adaptive filter for use in the echo cancellation system 
of Fig. 5. 



u DETAILED DESCRIPTION OF THE INVENTION 

s 

M- Referring now to Fig. 1 , a novel adaptive filter method is shown. A filter is chosen at 

□ 

m 110, preferably an adaptive finite impulse response (FIR) filter of an appropriate length N. Once 

the filter is chosen 110, convergence is achieved 120 using a convergence process, preferably a 

O 

(-* 20 least means square (LMS) approach. With convergence complete 120, the total energy of the 
filter coefficients is computed 130 and two energy thresholds, a first energy threshold and a 
second energy threshold, are set 140 relative to the computed total energy. Using the set energy 
thresholds, a subset of the converged filter coefficients is identified 150 by relating the starting 
filter coefficient with the first energy threshold and the end filter coefficient with the second 

25 energy threshold. A first filter delay, in accordance with the determined start coefficient, is 

incorporated 160 into the filter operation. In the filtration process, the converged subset of filter 
coefficients are then applied 170 to the input signal to create a modified signal used to filter 180 
the data through a summation or subtractor circuit in communication with the signal to be 
filtered. The error between the modified signal and signal to be filtered is monitored 190 and 

30 used to adapt 195 the filter coefficients and the process is repeated beginning with step 120. As 
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shown in Fig. 2, the adaptive filter system 200 comprises a filter 210 having a filter input 212 for 
receiving inputs from a far-end source 205, a filter output 215 for outputting a filtered result to a 
summation device 225, and an error input 230 for receiving an error signal generated by the 
comparison of an estimated signal against the actual signal. 
5 The present system will be described in greater detail by specific reference to a finite 

impulse response (FIR) filter. An FIR system is generally described by the difference equation: 

N 

y k = 2 hix k _i 

i = 0 

10 where y k represents an output signal that is the summation of the convolution of hi with input 
signal x k -i. One of ordinary skill in the art would appreciate that there are several methods for 
implementing an FIR system, including direct form, cascade-form, frequency-sampling, and 
lattice realizations. While the embodiments of the present invention will be described by 
reference to a direct form realization, the invention is not limited to direct form realizations and 



* 1 5 encompass any realization that could be effectively utilized in accordance with the teachings 
ry provided herein. Furthermore, while FIR filters are preferred for use in the present invention, a 

FIR filter form is not necessary and other filter forms, such as infinite impulse response (IIR) 

M ; filters, could be used. 

B 

las Designing an FIR filter requires the determination of N coefficients from a specification 

20 of the desired frequency response of the FIR filter. One of ordinary skill in the art would 
appreciate how to design a FIR system for a particular application and, more specifically, to 
define the appropriate length of the FIR filter. In the case of FIR filters being designed for echo 
cancellation applications, the desired frequency response of the FIR filter is to mimic the 
behavior of a cross-coupling transformation of a far-end signal, namely the nature and extent of 
25 the impulse response generated by the echo path. As such, the length of the FIR filter, 

designated by N, is calculated based upon the length of the impulse response that needs to be 
cancelled. For line echo cancellation, the impulse response is typically specified to be 64 ms in 
duration. The FIR filter length can therefore be approximated by multiplying the expected 
duration of the echo response, i.e., 64 ms, by the standard sampling rate in telephony systems, 8 
30 kHz. In this case, the FIR filter would comprise 512 coefficients or 512 taps. 

The taps, also referred to as filter coefficients or coefficients, of the FIR filter can be 
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realized by using a convergence process that employs any one of several known approaches. To 
achieve convergence, the approach would typically have to derive a set of coefficients that, when 
applied to an input signal, would generate a signal that was sufficiently close to (converged with) 
the signal that needed to be filtered. In the aforementioned echo cancellation application, the 
5 approach would be to update the coefficients once every sample, or every 125 microseconds. 
The more coefficients that need to be updated, the more computational resources and memory 
that are needed by the filter system. Reducing the number of coefficients would therefore reduce 
the computational complexity and memory requirements of the adaptive filter. This reduction 
becomes even more critical when dealing with applications requiring numerous adaptive filter 
10 operations, such as multi-channel echo cancellation devices. 

A conventionally used approach to achieve convergence comprises the LMS approach. 
One of ordinary skill in the art would appreciate the calculation that needs to be conducted in 
y order to conduct a convergence process using the LMS approach. Conceptually, by effectuating 
60 convergence, the LMS approach enables the determination of values for a set of coefficients that 

w 

py 15 comprise a transfer function that best approximates the transformation of an input signal to the 
u echo-containing output signal via the cross-coupling pathway. Depending on the length of the 

M= FIR filter, the LMS approach can be computationally complex and therefore comprises the most 

£3 

y, computationally intensive part of any filtration process. Additionally, the LMS approach 

5 exhibits slow convergence and requires the trial-and-error determination of an adaptation 

O 

M= 20 coefficient, denoted by |i, which controls the speed of convergence and, if improperly selected, 
effects the stability of the convergence calculation. 

Other approaches can be used and may be selected based upon a balance of factors 
including the rate of convergence, computational requirements, stability, and other properties. 
The rate of convergence is defined as the number of iterations required for the convergence 
25 process, in response to a set of inputs, to converge to the optimum solution. The computational 
requirements of the convergence process includes: (a) the number of operations (i.e., s 
multiplications, divisions, additions, and subtractions) required to perform one complete iteration 
of the convergence process; (b) the amount of memory needed to store the convergence process 
and accompanying data; and (c) the engineering investment required to program the convergence 
30 process. 
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As with other recursive algorithmic approaches, the LMS approach computes the 
coefficient values from the error signal by starting from some set of initial conditions that 
preferably approximates an assumed initial state of the system and iteratively configuring the 
coefficients to minimize the error signal when calculated in some mean-squared sense. 

More specifically, given an input vector u , and a vector of filter coefficients or 
weights w , then the minimum mean-squared error function can be written as: 

J(w) = a] -w H p-p H w+Rw, 
where a 2 d is the variance of the desired signal d(k), R is the auto-correlation matrix of the input 
signal u , and p is the cross-correlation between the desired signal d(k) and input u . 
The minimum value of J(w ) is: 

mmJ{w)~(j] -p H Rp 

for 

Wopftoio/ =R P- 

Here, w opljmal are the optimal weight of the filter in mean-squared sense. In LMS, the following 
update is used for the filter weights: 

0) 

where VJ is the gradient of J: 

VJ = 2E[u k " {d(k)-w[u k }] 
where E denotes a statistical expectation. The term in braces in the above equation is the error 
between the desired and estimated signal, which can be defined as: 

e(k) = {d(k)-w[u k ) 

In LMS, the statistical expectation is estimated by the instantaneous value of the gradient. 
Therefore: 

=w k + /je(ri)u k , 

(2) 
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It has been shown that the proper choice of ja should be: 

2 

0< M <- — , 

^max 

where A max is the maximum eigenvalue of the auto-correlation matrix R. Since R is not 
known and, therefore, A max is not known, one cannot necessarily choose a good value of \x. In 
5 practice, a value for \i is usually chosen by trial-and-error. The value of \x affects the filter 
performance. Smaller values of |i give higher signal-to-noise ratio but take more time to 
converge. Usually, a designer starts with a relatively large value of p, for fast initial 
convergence, and then chooses a smaller value for high SNR. 

In a preferred embodiment, after convergence is achieved on the FIR filter using the 
M- 10 above-described LMS method, the total coefficient energy is then computed. Total energy can 
□ be calculated based upon the coefficients, represented by w as calculated above, in accordance 



E = E Wi 2 



with the following equation: 

N-1 
0 

j\ where the energy calculation is performed over a set of i samples. As is known by those 

p skilled in the art, the total energy of a signal can be computed through other equations, such as: 

*J3 N-1 



i = 0 

The present invention is not limited to a single method for calculating energy and can be 
effectuated through any total energy calculation approach. 

In certain applications, such as echo cancellation, it is observed that energy levels vary 
over the course of a signal to be filtered and, as between signals, vary in profile and distribution. 

25 Because filter systems need to be designed to be long enough to filter the most substantial 

undesired frequencies likely to be received, they are often longer than necessary to account for 
many of the actual undesired signals encountered. More specifically, when dealing with line 
echo cancellation, it is known that an effective echo cancellation system should account for, at 
minimum, 32 milliseconds (ms) of delay. It is observed, however, that a substantial number of 

30 hybrid impulse responses do not have significant energy levels extending throughout the typical 
5-15 ms duration and, therefore, where sampling and coefficient calculation has occurred to 
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obtain coverage of a 5-15 ms duration echo, a majority of energy resides in only a subset of the 
total impulse response samples taken. 

As an example, a typical hybrid impulse response, caused by a specific far-end signal and 
sampled at 8000 Hz, may have a majority of the signal energy concentrated in the first few or the 
5 last few samples. A conventional filter having a length equal to 512 coefficients would therefore 
have a majority of its coefficients corresponding to energy levels at or near zero. Consequently, 
the filter process would use memory and computational resources for calculating and storing 
values that does not significantly improve system performance, reduce the signal to noise ratio, 
or improve filter quality. While possible to simply fix the filter at a smaller length, a smaller 
10 fixed filter length would prevent the system from being able to effectively respond to changes in 
the impulse response. 

P Referring back to Fig. 1, once convergence is achieved 120 and the total energy is 

calculated, two energy thresholds, a first energy threshold and a second energy threshold, ti and 
CO th respectively, are set 140 relative to the computed total energy. The values for t| and th denote 

ffj 1 5 the boundary that defines the selected portion of total energy to be filtered. The thresholds are 
~- calculated by applying a lower-end, pi, and higher-end, ph, cutoff, in the form of percentages, to 

E 

M» the total energy E. In echo cancellation applications, it is preferred that substantially all of the 

o 

u impulse response energy be filtered (approximately 99.99% of E) and, as such, the values for pi 

jO and ph are 0.001% and 99.991%. One of ordinary skill in the art would appreciate, however, that 

tax? 

jM- 20 the present invention is not limited to specific values for pi or ph. Those values can be selected 
by balancing quality considerations (requiring a larger portion of E to be filtered and therefore 
for pi «< ph) with the computational resources available (requiring fewer coefficients to be 
considered and therefore the difference between pi and ph need not be as great). The values for ti 
and th are then calculated in accordance with the following equations: 
25 ti = pi*E 

th = Ph * E 

Therefore, the two energy thresholds are approximately 0.001% of total energy E and 
approximately 99.99% of total energy E. The two thresholds bound 99.99% of the total energy 
to be filtered. 

30 With the two energy thresholds calculated, the filter of length N can be truncated in 
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operation to only address a subset of all N coefficients corresponding to the energy bounded by t| 
and th. Using ti and t h , corresponding start and an end coefficients are calculated by the 
following algorithm: 

Initial conditions: 
5 E t = 0 and k = 1 

Calculation: 

While E t < ti 

k-1 

E t =E t +E Wi 2 
10 i = 0 

k = k+ 1 

endwhile 

Start Coefficient = k 
P While Et + ti <t h 

W 

63 15 k-i 
y Et^Et+EWi 2 

O k = k+l 

y> endwhile 

n 

f 20 End Coefficient = k 

y3 Start_Coefficient and End Coefficient denote the coefficient boundaries associated with 

O 

y, the calculated energy boundaries. 

A filter delay, in accordance with the determined start filter coefficient, is incorporated 
into the filter operation. The delay can be incorporated through any method or system capable of 

25 accurately bypassing coefficients approximately less than the Start_Coefficent and initiating 
coefficient usage around or at Start_Coefficient. An exemplary method would comprise 
adjusting a filter data pointer to a position associated with Start_Coefficient or slightly less than 
Start_Coefficient. In the course of operation, after bypassing initial coefficients to reach at or 
about Start_Coefficient, the filter applies the coefficients approximate to or greater than 

30 Start_Coefficient and less than or approximate to End_Coefficient to model the signal to be 

filtered, terminating at or around End_Coefficient. Where pi «< ph, the application of a subset 
of total coefficients delivers substantially similar filtration quality to the application of all 
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coefficients, but with fewer computational and memory resources required. 

In operation, the coefficients are converged, total energy is calculated, the lower and 
upper energy thresholds are calculated, the corresponding start and end coefficients are 
determined, and, with those start and end coefficients, a subset of all coefficients is identified. 
5 As shown in Fig. 1, the identified subset is then used to filter a signal by applying the identified 
subset of coefficients to an input signal, generating a modified input signal, and applying that 
input signal to the signal to be filtered in order to remove the undesired frequencies modeled by 
the application of the coefficient subset to the input signal. One of ordinary skill in the art would 
appreciate that the filter process can be conducted in numerous ways and that the present 
1 0 invention is not limited to a specific method for using the identified coefficient subset for 
filtration. 

Referring now to Fig. 3, a second embodiment of the present invention provides for 
selection 310 of a filter, preferably an adaptive finite impulse response (FIR) filter of an 
appropriate length N. Once the filter is chosen 310, convergence is achieved 320 using a 
pi 15 convergence process, preferably a least means square (LMS) approach. With convergence 
O complete 320, the process to identify a subset of converged coefficients is initiated 325. A set of 

5 

\a coefficients (Selected Coefficients), which represents a sub-set of the entire set of converged 

H coefficients and is preferably pre-designated but, optionally, may be dynamically chosen, is used 

to calculate 330 the average energy of the portion of the signal corresponding to the Selected 
I 20 Coefficients. The average energy is then scaled using a scaling factor to produce 340 a threshold 
(herein referred to as a threshold value) against which the rest of the coefficient are to be 
evaluated. 

Starting from a specific coefficient, preferably the first coefficient, the average energy 
value of a series of coefficients, where the series number is equal to the number of Selected 

25 Coefficients, is calculated 350. This value is then compared 360 against the previously 

calculated threshold value. If the previously calculated threshold value is greater than 363 the 
value calculated by step 350, the next coefficient series, once again where the series number is 
equal to the number of Selected Coefficients, is selected 367 and evaluated in the same manner 
as steps 350 and 360. If the previously calculated threshold value is equal to or less than 369 the 

30 value calculated by step 350, the first coefficients of the subset of converged coefficients have 
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been identified 370. In parallel or in series, the last coefficients of the subset of converged 
coefficients can be identified 380 in a similar manner by performing the same calculation, except 
beginning with the end coefficient and having the series move toward the first coefficient. 
Alternatively, the end coefficient can be identified by adding a fixed number of coefficients to 
5 the start coefficient to yield an end coefficient. 

The selection of the Selected Coefficients can be performed by any method, most 
preferably using those methods that are appropriate to the application to which the adaptive filter 
is being applied. For example, the Selected Coefficients can be a pre-designated set of 
coefficients with fixed start and end coefficients that may correspond to a reliably present energy 
10 presence in the signal to be filtered. Alternatively, the Selected Coefficients can comprise a 

fixed start coefficient with variable end coefficients depending on the signal strength of the fixed 
start coefficient. Where computational resources and time permit, the Selected Coefficients can 
also be dynamically chosen by selecting the set of coefficients having the highest signal strength. 
One of ordinary skill in the art would appreciate that nearly any method can be used to identify 
15 the Select Coefficients and such methods are preferably unique to the adaptive filter application. 
O More specifically, the Selected Coefficients are preferably selected based upon a best estimate as 

s 

M= to what area, in the signals to be filtered, would typically comprise noise. While signals can vary 

D 

Tl greatly, within a particular application area, signals often have at least one area that could be 

predictably assumed to predominantly comprise noise. For example, within echo cancellation 
jjt 20 applications, the end coefficients of a typical line echo tend to predominantly comprise noise due 
to propagation delay. Therefore, it is preferred to select a portion of such end coefficients as the 
Selected Coefficients. It is further preferred to use between 10 and 20 (more preferably 16) of 
the last end coefficients. One of ordinary skill in the art would appreciate that which coefficients 
are selected, and the number of coefficients selected, depends upon the nature of the application 
25 and the type of signals expected. 

With both the start and end coefficients identified 370, 380, a first filter delay, in 
accordance with the determined start coefficient, is incorporated 390 into the filter operation. In 
the filtration process, the converged subset of filter coefficients are then applied 392 to the input 
signal to create a modified signal used to filter 394 the data through a summation or subtractor 
30 circuit in communication with the signal to be filtered. The error between the modified signal 
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and signal to be filtered is monitored 396 and used to adapt 398 the filter coefficients and the 
process, beginning with step 320, is repeated. 

In a preferred embodiment, after a filter is chosen, preferably a FIR filter, and 
convergence is achieved, preferably using the LMS method, the average energy of the Selected 
Coefficients is calculated. The average coefficient energy for the Selected Coefficients series is 
computed. Average energy can be calculated based upon the Selected Coefficients, represented 
by w as calculated above, in accordance with the following equation: 

M 

Average E = I Wj 2 

i = X 

M-x 

where the energy calculation is performed over a set of i samples. As previously stated, the 
average energy of a signal can be computed through other equations and the present invention is 
O not limited to a single method for calculating average energy. Where the last 16 taps of a filter 
©15 of length N=5 12 are selected, as is preferred for echo cancellation applications, the values for x 

and M are 496 and 512 respectively. 
O With the average energy calculated for the Selected Coefficients, the average energy is 

s 

\a then scaled using a scaling factor. The scaling factor is experimentally determined in accordance 

Q 

jj, with the nature of the application for which the adaptive filter is being used. Higher scaling 

€? 20 factors generate higher threshold values, thereby providing for less echo cancellation, lower filter 

O 

LI quality, and higher computational and memory resource savings. Lower scaling factors generate 

lower threshold values, thereby providing for higher echo cancellation, greater filter quality, and 
lower computational and memory resource savings. One of ordinary skill would appreciate the 
scaling factor is dependent upon the type of application to be filtered. In echo cancellation 
25 applications, it is preferred to use a scaling factor of 4. The threshold value (TV) is therefore 
average energy (AvgE), as calculated above, times the scaling factor (SF): 

TV = AvgE * SF 

With TV calculated, each series of coefficients, where a series is equal to M-x (the 
number of Selected Coefficients), is evaluated by calculating the series' average energy times the 
30 scaling factor. While one can begin the series at any place, particularly if it is known that a 
signal tends to appear at or around specific coefficients, for echo cancellation applications, it is 
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preferred to begin with the first coefficient and move up an incremental M-x number of 
coefficients in each series. By selecting windows of coefficients, one can incrementally compare 
each series against the threshold value to determine whether a portion of the converged 
coefficient set exceeds a base noise energy level and, therefore, represents a signal to be filtered. 
If so, the coefficients corresponding to that set are relevant to an on-going filtration process. If 
not, the coefficients are not relevant and are using computational and memory resources 
unnecessarily. By basing the comparison process against a known signal area to contain noise, 
one can have a signal to noise ratio based method of truncating a filter to efficiently filter an 
incoming signal. 

Assuming that the beginning point is the first coefficient, k is equal to the number of 
Selected Coefficients, SC is equal to a start coefficient, and TV is equal to a previously 
calculated threshold value: 

Initialization 

SAvgE - 0; SC=0 

Calculation: 

While TV > SAvgE 

SAvgE = 0;SC = SC + k 

SC + k 

SAvgE 2 Wi 2 

i = SC 

k 

endwhile 

Start_Coefficient - SC 

Start_Coefficient is the first coefficient in the tested series. As an example, assuming the 
scaled average energy of coefficient series 18 to 34 exceeded the threshold value, the start 
coefficient (SC or Start_Coefficient) would be 18. One of ordinary skill in the art would 
appreciate that protection mechanisms could be incorporated into the above calculation to protect 
against situations where, for example, the requirement TV < SAvgE is not met and SC begins to 
exceed the total length of the filter. 

The end coefficient can be identified in the same manner, except beginning with the last 
coefficient and moving the series backward by a number of coefficients equal to the pre-selected 
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number of coefficients. For example, the above calculation would begin with SC = 5 12 
(assuming a filter of length N=512) and then incrementally moving downward by 16 (assuming 
k = 16) through a modification of SC = SC + k to SC = SC -k. Alternatively, the end coefficient 
could be a fixed number of coefficients after the start coefficients, assuming the fixed number 
5 yields an end coefficient less than or equal to the total filter length. 

With the start coefficient and end coefficients identified, the coefficient boundaries are 
established and the subset of coefficients for use in the filtration process have been determined. 
A filter delay, in accordance with the determined start coefficient, is incorporated into the filter 
operation. The delay can be incorporated through any method or system capable of accurately 
10 bypassing coefficients approximately less than the start coefficient and initiating coefficient 

usage around or at the start coefficient. An exemplary method would comprise adjusting a filter 
data pointer to a position associated with the start coefficient or slightly less than the start 
coefficient. In the course of operation, after bypassing initial coefficients to reach at or about the 
start coefficient, the filter applies the coefficients approximate to or greater than the start 
JJJ 1 5 coefficient and less than or approximate to the end coefficient to model the signal to be filtered, 
O terminating at or around the end coefficient. 

y, In operation, the coefficients are converged, energy is calculated for the Selected 

P Coefficients, the energy threshold is calculated, the corresponding start and end coefficients are 

yj determined, and, with those start and end coefficients, a subset of all coefficients is identified. 

n 

«» 2Q The identified subset is then used to filter a signal by applying the identified subset of 

coefficients to an input signal, generating a modified input signal, and applying that input signal 
to the signal to be filtered in order to remove the undesired frequencies modeled by the 
application of the coefficient subset to the input signal. One of ordinary skill in the art would 
appreciate that the filter process can be conducted in numerous ways and that the present 
25 invention is not limited to a specific method for using the identified coefficient subset for 
filtration. 

The above-described novel adaptive filter methods and systems can be effectively 
deployed in a telecommunications system in the form of novel echo cancellation methods and 
systems to effectuate high quality communications, particularly as between users of a public 
30 switched telephone network (PSTN) and users of a packet based network (e.g., the Internet). 
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Referring to Figure 4, a telecommunication system 400 is shown comprising a PSTN 410 having 
a plurality of telephonic systems 415 (e.g., telephones and fax machines) and a packet based 
network 420 having a plurality of networked systems (e.g., file servers 430, email servers 440, 
computers 450) linked via routers 460. Mediating between the two networks is a gateway 470 
comprising a plurality of echo cancellation devices capable of substantially reducing echo 
generated by the movement of received inputs through various cross-coupling pathways, and a 
plurality of digital to analog and analog to digital encoders and decoders . Because of the 
potential volume of signals traveling through such a system 400, echo cancellation can be 
occurring on numerous channels, carrying different signals, concurrently. It is therefore 
important to optimize echo cancellation in a manner that reduces the computational complexity 
and memory requirements associated with the echo cancellation process. The novel adaptive 
filter method and system provided herein can be used to achieve such a result. 

Referring to Fig. 5, a far-end signal x 510 from a far-end source 515 is received locally at 
a communication input 520. As a result of the previously noted imperfections in the local 
system, a portion of the signal x 510 is echoed back to a transmit line 530 via a cross-coupling 
path 525 that intersects with the transmit line 530 at an input 535. In the course of engaging the 
cross-coupling path 525, whether an acoustic or line echo pathway, the signal x 510 undergoes a 
transformation, as a function of an impulse response, that modifies the original signal x 510 into 
a new signal s 540, referred to herein as the echo response 540. The echo response 540 is 
illustrated here as a signal s 540 corresponding to the following equation: 

s=h*x 

where h is the impulse response of the echo characteristics. 

Also being communicated through the transmit line 530 is a near-end signal v 545, 
communicated from a near-end source 550. The input signal v 545 combines with the echo 
response s 540 to generate a combined signal y 555. Therefore, the signal sent from the near-end 
source 550 to the far-end receiver 580, absent echo cancellation, is the signal y 555, which is the 
sum of the near-end signal v 545 and the echo response s 540. 

To reduce and/or eliminate the echo response component s 540 from the signal y 555, a 
typical system uses an echo canceller 560 having a filter 565 that is capable of applying an 
impulse response fi, which is an estimate of the actual impulse echo response h experienced by 
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the far-end signal x 510 as it engages the cross-coupling path 525. As such, a further signal s 
570 representing an estimate of echo response s 540 is generated by the echo canceller 560 in 
accordance with the following equation: 

s = h *x 

5 The echo canceller 560 subtracts the echo estimate signal s 570 from the signal y 555 to 

generate a signal e 575 that is returned to the far-end receiver 580. The signal e 575 thus 
corresponds to the following equation: 

e = s+v- s » v 

The signal returned to the far end receiver 580 is therefore dominated by the signal v of 
10 the near-end source 550. To the extent the impulse response h more closely correlates to the 
actual echo impulse response h, then s 570 more closely approximates s 540, resulting in the 
minimization of the magnitude of the echo signal component s 540 on the signal e. 

C 

O An adaptive filter 565 is used to generate the echo signal component s 570. In its 

yy 

J simplest form, the adaptive filter 565 generates an echo estimate, i.e., s 570, by obtaining 

W 1 5 individual samples of the far-end signal x 5 1 0 on a receive path 513, convolving the samples 

fy 

Q with an impulse response model of the system, i.e., h, and then subtracting, at the appropriate 

^ time, the resulting echo estimate, s 570, from the received signal y 555 on the transmit channel 

Q 530. The conventional adaptive filter is a FIR filter using a LMS method for achieving tap 

convergence. 

^ 20 The novel adaptive filter methods and systems presented herein can be used to improve 

the calculation of the echo impulse response h by, among other things, reducing the 
computational complexity and memory requirements of the coefficient calculation conducted 
within the adaptive filter. Shown in Fig. 6, an embodiment of the novel filter 665 of the present 
invention is used to generate the echo signal component s 670. After having achieved 
25 convergence on a FIR filter and truncated the filter, in accordance with either one of the 

previously described methodologies, the adaptive filter 665 generates an echo estimate, i.e., s 
670, by obtaining individual samples of the far-end signal x 610 on a receive path 613, 
convolving the samples with the calculated coefficients, and then subtracting, at the appropriate 
time, the resulting echo estimate, § 670, from the received signal y 655 on the transmit channel 
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630. Ongoing adaptation of the filter occurs by the adjustment of the filter coefficients, 
represented by the arrow 690 extending through element 665. 

To provide a delay to insure that only coefficients approximate to or greater than the SC 
are employed, a delay is applied where D - 1 is a specific value of delay related to the value of 
5 SC. The start and end coefficients may have been calculated using either of the first or the 

second novel adaptive filter methods and systems described above. In a preferred embodiment, 
D is equal to SC so the delay is equal to SC - 1 . One of ordinary skill in the art would appreciate 
that different values of a delay, as related to SC could be used and depends upon how much 
adherence one of ordinary skill in the art would want the designed system to have relative to the 
10 calculated SC. Some may prefer to strictly adhere to the calculated threshold (delay = SC -1). 
Others may prefer to incorporate a greater buffer (delay « SC) in order to insure greater quality 

H 3 while sacrificing some of the potential computational and memory resource savings. 

O 

O The signal s 670 is produced as a function of the transfer function denoted by z" D 685. 

8 a I 

P The values of D extend from the point after the initial delay up to a point proximate to the EC. 

W 15 In a case where the delay is equal to SC - 1 and the termination point is chosen to equal EC, 

RJ 

□ keeping in mind the termination point could be slightly greater than or less than EC, depending 

jj\ upon design considerations, the transfer function z" D 685 can be modeled as: 

O z" sc Z ' (sc+1) z" (sc + 2) z~ (sc + 3) ...z~ (EC ~ 3) z' (EC ~ 2) z (EC ~ l) z" EC 

fc j3 Operationally, the novel echo cancellation application has achieved superior performance 

!*f 20 results. Applying an echo cancellation system using the first novel adaptive filter to a 

conventional hybrid response, a FIR filter of length N=512 was chosen and converged according 
to the standard LMS approach. Using this conventional approach, the signal to noise ratio was 
48 dBs. Using the approach of the present invention, a FIR filter of length n=5 12 was chosen, 
converged using a LMS approach, and subjected to a low and high threshold, set at 0.001% and 
25 99.99% of the total energy. The resulting threshold application bounded the start and end 
coefficients at 51 and 127 respectively, reducing the filter length to 77. The filter result was 
acceptable, producing a signal to noise ratio of 42 dBs. Without an unacceptable decrease in 
performance, the present invention provided for a substantially reduced filter length and 
therefore a reduced computational and memory resource burden. 
30 The present adaptive filter method and system can be employed in numerous applications 
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employing adaptive processes in conjunction with convolutional coding. Accordingly, another 
embodiment of the present invention includes a novel method and system for channel 
equalization. Equalizers are a class of communication system devices used to compensate for 
distortion experienced in communication channels. Fixed equalizers have the average electrical 
characteristics of the channel pre-determined and a fixed amount of equalization is therefore 
designed into the equalizer to compensate for the distortion. Adjustable equalizers monitor the 
channel and provide for equalization that varies, as necessary, to match the distortion determined 
at the time of monitoring. Adjustable equalizers, also known as adaptive equalizers, can provide 
for the adaptation using an adaptive process to generate a modified signal, capable of canceling 
out distortions in the signal, or an equalizing signal. 

Applied to a channel equalization application, the present invention is used to generate 
the equalizing signal by achieving convergence on a FIR filter, truncating the filter in accordance 
with either the first or second novel adaptive process disclosed herein, and generating a 
equalizing signal by obtaining samples of a received signal and convolving the samples with the 
calculated coefficients. The channel equalizer then applies, at the appropriate time, the resulting 
equalizing signal to the received signal. Ongoing adaptation of the filter occurs by utilizing an 
error signal, generated from a comparison of the equalizing signal to the received signal, and 
adjusting the coefficients based upon the error signal to obtain convergence. 

To provide a delay to insure that only coefficients approximate to or greater than the start 
coefficient (hereafter abbreviated as SC) are employed, a delay is applied where D - 1 is a 
specific value of delay related to the value of SC. The start and end coefficients may have been 
calculated using either of the first or the second novel adaptive filter methods and systems 
described above. In a preferred embodiment, D is equal to SC so the delay is equal to SC - 1 . 
One of ordinary skill in the art would appreciate that different values of a delay, as related to SC 
could be used and depends upon on how much adherence one of ordinary skill in the art would 
want the designed system to have relative to the calculated SC. Some may prefer to strictly 
adhere to the calculated threshold (delay = SC -1). Others may prefer to incorporate a greater 
buffer (delay « SC) in order to insure greater quality while sacrificing some of the potential 
computational and memory resource savings. 
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The equalizing signal is produced as a function of the transfer function denoted by z~ D . 
The values of D extend from the point after the initial delay up to a point proximate to the end 
coefficient (hereafter abbreviated as EC). In a case where the delay is equal to SC - 1 and the 
termination point is chosen to equal EC, keeping in mind the termination point could be slightly 
greater than or less than EC, depending upon design considerations, the transfer function z* D 785 
can be modeled as: 

Z -SC Z -(SC+1) z -(SC + 2) z -(SC + 3) z -(EC-3) z <EC-2) Z -(EC~1) ^EC 

The present methods and systems provide for an adaptive filter that significantly reduces 
the memory requirement, memory bandwidth, and computational resources necessary to operate 
the filter. Applied to the problem of echo cancellation, embodiments of the present invention 
obtained a reduction of computational resource usage while maintaining an acceptable 
performance level. While various embodiments of the present invention have been shown and 
described, it would be apparent to those skilled in the art that many modifications are possible 
without departing from the inventive concept disclosed herein. For example, the adaptive filter 
has been applied in an echo cancellation application. It would be appreciated by one of ordinary 
skill in the art that the filter can be used in any application where the convergence, stability, 
computational requirements, and memory bandwidth characteristics of the novel filtration 
method and system could be effectively applied. 
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